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The knowledge of the redshift of multiple images in cluster-lenses allows to determine precisely 
the total projected mass within the Einstein radius. The observation of various multiple im- 
ages in a same cluster is opening new possibilities to constrain the curvature of the universe. 
Indeed, although the influence of Qm and JIa on the images formation is of the second order, 
observations of many multiple images at different redshifts formed by a regular cluster-lens 
should allow to constrain very accurately the mass distribution of the cluster and to start to be 
sensitive to the cosmological parameters entering the diameter angular distances. We present, 
analytical expressions and numerical simulations that allow us to compute the expected error 
bars on the cosmological parameters provided an HST/WFPC2 resolution image and spectro- 
scopic redshifts for the multiple images. Numerical tests on simulated data confirm the rather 
small uncertainties we could obtain this way for the two popular cosmological world models: 
= 0.3±0.24, D.X = 0.7±0.5 or fi™^ = l.±0.33, D.x = 0.±1.2. 



1 Introduction 

Recent works on constraining the cosmological parameters using the CMB and the high redshift 
supernovae seem to converge to a "standard cosmological model" favouring a flat universe with 

~ 0.3 and ^Ix ~ 0.7: White!. However these results are still uncertain and depend on some 
physical assumptions, so the flat Om, = 1 model is still possible (Le Dour et al.%). It is therefore 
important to explore other independent techniques to constrain these cosmological parameters. 

In cluster gravitational lensing, the existence of multiple images - with known redshifts - 
given by the same source allows to calibrate in an absolute way the total cluster mass deduced 
from the lens model. The great improvement in the mass modeling of cluster-lenses that includes 
the cluster galaxies halos (Kneib et alM, Natarajan & KneibH) leads to the hope that clusters 
can also be used to constrain the geometry of the Universe, through the ratio of angular size 
distances, which only depends on the redshifts of the lens and the sources, and on the cosmolog- 
ical parameters. The observations of cluster-lenses containing large number of multiple images 



lead Link & PierceQ (hereafter LP98) to investigate this expectation. They considered a simple 
cluster potential and on-axis sources, so that images appear as Einstein rings. The ratio of such 
rings is then independent of the cluster potential and depends only on and assuming 
known redshifts for the sources. According to them, this would allow marginal discrimination 
between extreme cosmological cases. But real gravitational lens systems are more complex con- 
cerning not only the potential but also off-axis positions of sources. They conclude that this 
method is ill-suited for application to real systems. 

We have re-analyzed this problem building up on the modeling technique developed by us. 
As demonstrated below, we reach a rather different conclusion showing that it is possible to 
constrain and Q\ using the positions of multiple images at different redshifts and some 
physically motivated lens models. 



2 Influence of and Q\ on the images formation 

2. 1 Angular size distances ratio term 

In the lens equation: Os = Oi — pr^V(^^(6'i), the dependence on and is solely con- 

^os 

tained in the term F = DqlDls I Dqs- For a given lens plane, F{zs) increases rapidly up to 
a certain redshift and then stalls, with significant differences for various values of the cosmo- 
logical parameters (see Fig. [l|). Thus in order to constrain the actual shape of F{zs) several 
families of multiple images are needed, ideally with their redshifts regularly distributed in F[zs) 
to maximize the range in the F variation. 
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Figure 1; Left. F{zs) for zi = 0.3 and various cosmological models. Right. Fiza'i) I F(za\) as a function of f2m 

and JIa for zi = 0.3, z^i — 0.7 and Zs2 = 2. 

If we consider fixed redshifts for both the lens and the sources, at least 2 multiple images are 
needed to derive cosmological constraints. In that case F has only an influence on the modulus 
of 6i — 6s- So taking the ratio of two different F terms provides the intrinsic dependence 
on cosmological scenarios, independently of Hq. A typical configuration leads to the Fig. |l| 
plot. The discrepancy between the different cosmological parameters is not very large, less than 
3% between an EdS model and a flat low matter density one. The figure also illustrates the 
expected degeneracy of the method, also confirmed by weak lensing analyzes, with a continuous 
distribution of background sources {e.g. Lombardi & BertinH ). 



2.2 Relative influence of the different parameters 



We now look at the relative influence of the different parameters, including the lens parameters, 
to derive expected error bars on Q,^ and To model the potential we choose the mass density 
distribution proposed by Hjorth &: KneibS, characterized by a core radius, a, and a cut-off radius 
s 3> a. We can then get the expression of the deviation angle modulus Dq^ =|| Oj — Os ||. 

For 2 families of multiple images, the relevant quantity becomes the ratio of 2 deviation 
angles for 2 images On and 6j2 belonging to 2 different families at redshifts Zgi and Zs2- Let's 

define Roji = — With several families, the problem is highly constrained because a single 

^012 

potential must reproduce the whole set of images. In practice we calculate ^i^^^i^ versus the 
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different parameters it depends on. For a typical configuration: zi = 0.3, Zgi = 0.7, Zs2 = 2, 

- — = 2, — = 10 {6a = a/Doifis = s/Dol) and we assume Vim = 0.3 and Q.\ = 0.7. We then 
"n "a 

obtain the following orders of magnitudes for the different contributions : 
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As expected, even with 2 families of multiple images the infiuence of the cosmological param- 
eters is of the second order. The precise value of the redshifts is quite fundamental, therefore 
a spectroscopic determination {dz = 0.001) is essential. The position of the (fiux-weighted) 
centers of the images are also important. With HST observations we assume dOj = 0.1". 

So even if the problem is less dependent on the core and cut-off radii, they will represent 
the main sources of error. Taking dOa/Oa = dOg/Os = 20 %, we then derive the errors d^lrn and 
d^lx from the above relation in these two cosmological scenarios : 

J)m = 0.3 ±0.24 J]a = 0.7±0.5 or 0^ = 1 + 0.33 fix = ±1.2 

As confirmed by the Fig. [l| degeneracy plot, the method is more sensitive to matter density 
than to the cosmological constant. 



3 Constraint on {Qm,^x) from strong lensing 

3.1 Method and algorithm for numerical simulations 

We consider basically the potential introduced in section 2.2. After considering the lens equation, 
fixing arbitrary values (0,^,0,^) and a cluster lens redshift z/, our code can determine the images 
of a source galaxy at a redshift Zs- Then taking as single observables these sets of images as 
well as the different redshifts, we can recover some parameters (the more important ones being 
(TO) 6a or 6s) of the potential we left free for each point of a grid (0^,1^^). The likelihood of the 
result is obtained via a ^'^-minimization, where the is computed in the source plane. 



3.2 Numerical simulations in a typical configuration 

To recover the parameters of the potential ( ie cq, 6a, 6s and adjusted lens parameters), we 
generated 3 families of images with regularly distributed source redshifts. 

For starting values (r2^,$7^) = (0.3,0.7) we obtained the Fig. |2| confidence levels. The 
method puts forward a good constraint, better on Qm than on Qx, and the degeneracy is the 
expected one (Fig. 2). Concerning the free parameters, we also recovered in a rather good way 
the potential, the variations being Actq ~ 150 km/s, A6a ~ 3" and A^^ ~ 20". 

This is an "ideal" case, of course, because we tried to recover the same type of potential we 
used to generate the images, the morphology of the cluster being quite regular and the redshift 



-40 -20 20 40 2 




0.0 0-2 0-4 0.6 OS 1.0 



-40 -20 20 +0 

Figure 2: Left. Generation of images by a 2; = 0.3 cluster with ctq = 1400 km/s, 6a ~ 13.54" and 6s = 145.8". 
Close to their respective critic lines, we see 3 families of images at Zsi — 0.6, Z32 = 1 and Zs3 — 2. Right. Solid lines 
: X^i^rn, confidence levels obtained for this configuration. Generating arbitrary values: (57^, ^1) = (0.3, 0.7). 
Dashed lines: x'^{^rn,0,\) confidence levels obtained considering 10 clusters in this same configuration. 



range of the sources being wide enough to check each part of the F curve. Such simple approach 
can be appHed to regular clusters like MS2137-23, which shows at least 3 families of multiple 
images including a radial one. But the spectroscopic redshifts are still missing for the moment. 

4 Conclusions & prospects 

Following the work of LP98, we discussed a method to obtain informations on the cosmological 
parameters and 0,x while reconstructing the lens gravitational potential of clusters with 
multiple image systems at different redshifts. 

This technique gives degenerate constraints, and ^Ix being negatively correlated, with a 
better constraint of the matter density. With a single cluster in a typical lensing configuration 
we can expect the following error bars : = 0.3ib0.24, 0,x = 0.7ib0.5. To perform that, several 
general conditions must be fulfilled: a cluster with a rather regular morphology, "numerous" 
systems of multiple images, a good spatial resolution (HST) and spectroscopic precision for the 
different redshifts that should be also regularly distributed, from zi to high values - this requires 
deep spectroscopy on 8- 10m class telescopes due to the faintness of the multiple images. 

Combining the study of about 10 different clusters would tighten the error bars and lead 
to meaningful constraints. The dashed line confidence levels in the Fig. H are the result of a 
numerical simulation made with 10 identical clusters. We are encouraged by more and more 
known observations including systems with multiple sources and we plan to apply in a first time 
this technique to clusters like MS2137-23, MS0440+02, A370, AC114 and A1689. 
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